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SUMMARY 


An asymptotic description is derived for the interaction of a weak normal 
shock wave and a turbulent boundary layer along a plane wall. In the case stud- 
ied the nondimens ional friction velocity is small in comparison with the nondi- 
mensional shock strength, and the shock wave extends well into the boundary 
layer. Analytical results are described for the local pressure distribution 
and wall shear, and a criterion for incipient separation is proposed. A compar- 
ison of predicted pressures with available experimental data includes the effect 
of longitudinal wall curvature. 


INTRODUCTION 


In transonic flows a transition from supersonic to subsonic speeds typi- 
cally occurs through a shock wave which is normal to a solid boundary. Across 
the boundary layer along this surface, the upstream fluid velocity decreases 
from its value in the external flow to zero at the wall, and the strength of an 
incident shock wave must decrease to zero in the supersonic part of the boundary 
layer; at the wall there is no discontinuity in pressure. 

In the undisturbed turbulent boundary layer, viscous forces are important 
only in a very thin wall layer having thickness of the order of the local vis- 
cous length and velocity of the order of the friction velocity. In most of the 
boundary layer the mean velocity profile is nearly uniform, with a decrease from 
the external-flow value of the same order as the friction velocity. The flow 
changes caused by a weak incident normal shock wave will depend on the relative 
sizes of this velocity variation and the difference < .ween the external-flow 
Mach number and one, i.e., on the ratio of the none itisional friction velocity 
to the nondimens ional shock-wave strength. Asympt^.ic flow descriptions have 
been derived for large (ref. 1), moderate (refs. 2,3,4) and small (ref. 5) 
values of this ratio, and the oblique-shock problem has also been studied (ref. 6). 
In each case the local mean pressure gradient and fluid deceleration are large 
enough that the flow near the shock wave may be described as an inviscid rota- 
tional flow in most of the boundary layer. In a thinner sublayer, changes in 
the turbulent stresses do influence the changes in velocity, and here some model 
must be chosen (e.g., a mixing-length model) for describing these stresses. T’.e 
displacement effect of the flow changes in this sublayer is small enough that 
the dominant terms in the pressure are not affected. Thus the pressure is cal- 
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culated from the inviscid-flow equations and then is substituted 
layer equations for the calculation of changes in wall shear. 


into the sub- 


The present work is a continuation of the work of reference 5. Asymptotic 
solutions are obtained in the limit as the nondimensional friction velocity and 
shock-wave strength approach zero, such that the ratio of frictlS velocity to 
shock strength also approaches zero, as does the distance from the wall to^he 

Jould Jn by the bo ^dary-layer thickness. The pressures 

{fZLt I! * " Ce ? a 5 e ahown to be edified by higher-order and curvature 
wf facts, and the calculation of wall shear is carried out, lending to a pro- 
posed criterion for incipient separation. P 


SYMBOLS 


K 

L 

Me 

P»P e »^o 

Re,R<? 

U,U e 

U T 

X,x 

Y »y»y 

Y 

6 

e 

K 

n 


skin-friction coefficient 
local wall curvature, nondimensional with L -1 
length of boundary layer, up to shock wave 
Mach number in external flow ahead of shock 

local static pressure, and external-flow static and stagnation 

pressures ahead of shock, nondimensional with sonic pressure ahead 
or shock 

Reynolds number based on L and on external-flow quantities or sonic 
values 

local velocity and external-flow velocity, nondimensional with 
critical sound speed 

nondimensional friction velocity i u ^ — U ^ r 
2 ‘ 2 e c f 

X/L, (M g - 1) 1 2 X/6 ; X » coordinate along wall 

Y/L, Y/6 . (M g 2 - 1) 1 y/u T ; Y « coordinate normal to wall 
ratio of specific heats 

boundary-layer thickness ahead of shock, nondimensional with L 
U e - 1 

von Karman constant, taken equal to 0.41 

pressure-gradient parameter in velocity profile 

wall shear stress, nondimensional with undisturbed value 
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PRESSURE DISTRIBUTION 


-.f Veloc ^l ty defect la y er « y “ Y/6 = 0(1), where 6 = 0(u T ), and ahead 

kr-1 n W3Ve U ~ 1 + e + u x u 01<y). a s in figure 1. We take u 01 = k~ 1 i n y 

K . 0 J/LJ ““ ' y) for 0 < y < 1 and a 0x - 0 for y > 1 (a.g., ref* 7), Ohara 
. ?‘ 41 and for * ero Pressure gradient n is about 0.5. The sonic line is lo- 

ve v .h /,,;.? 1 ? (2 " ; “V” “T « « this value la sZil it tl 

y thin wall layer Y - 0(u T lRe i), and the velocit” V = 0(u T ) has the same 

f°™ aa °fl an 1 i " com P ressible flow having density and viscosity coefficient 

4 311 Values in the actual flow. For tt “1 Re” 1 « Y « 6 a mixing- 
ength model is used for the shear stress, and the velocity profile obtained is 

profiles asy “P totlcall y with both the defect-layer and the wall-laye^ 

6 u and m3 i ng c ° ndibions Provide a relationship among the quantities 

IcZllnt < ! UCh relat±onshi P I® found from the momentum integral 
equation, in the manner described for incompressible flow, e.g., in reference 7 

sirssLT" ln retained - For si - piicity the e„ t haip, e n oce '■ 


For Y - 0(6) immediately downstream of the shock wave, the differential 
equations and shock jump conditions (expanded for u T - 0, e - 0, and u /e + 0) 

- on 5 hat P ressure changes 0(u T ) occur over a distance x * (H g 2 - l)-l/2 T x /6 

- 0(1). The mean velocity profile has a term (1 + e)-l plus terms 0(u T ) con- 

shn^" 8 3 ?° W ? rota f iona l P art and an unknown irrotational part (fig. 2). The 

j , wave s located at x = 0(u T /e), and so the shock jump conditions are ex- 

?iof L th^a^ahl ri6S ab r fc V FOr a first a PP-x Jtion LapUce's equa- 

sources alone x = v’ V JV? 2U/? 3 quarter plane by distributing fluid 
a °?“? 7 8 7 °\ For (x2 y 2 ) 1/2 greater than about 2, it is found that the 

solution is nearly that due to a point source which represents the change in 
boundary-layer dispiacement thickness (as also noted in refs. 3 . c rec- 

"v “'„fl le l°7 U ;: ° £ hi8het ° td "' *««»«“« numerically at rlllil 

of the mean flS* the 6 rCpresenbs chan S es vovticity along a streamline 

Hnec „ ! , ? he accounts for the difference between streamlines and 

enters^iii th^irwtaH Calcula J lon of the rotational part of U, and also 
enters in the irrotational part through the boundary condition at x = 0 The 
wall pressure with these effects ncluded is 


P /P 
w e 


1 + 2 ye + y (2v-i) e 2 


+ (1 + (2y-l)e } & u (i + 21 h. I ± 3-18 I ± 1.5 n 2 “ u 0i 
TTT 2 1C 1 + n ' X L ~T 


(n) dn 


2 2 
x + n 


(l) 


The solution (1) is logarithmically infinite as x ■> 0, as must be expected since 

the nonHn! P rassure change is 0(e) rather than 0(u T ). Numerical solution of 
e non. inear transonic small-disturbance equations would be required to obtain 

5?™/° 3 re8i ° n near X = 0 which is ^ u T /e << i lor a lln 

with longitudinal curvature, the. pressure gradient in the iJviscid flow behind a 
normal shock wave has a logarithmic singularity described by 


A p = { x 6 In 6 /x 2 + y 2 - y 6 (tan -1 ^ - y) } 


( 2 ) 
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The boundary- layer interaction effect for x = 0(1) is not large enough to change 
this expression, and so a curvature term can simply be added as a correction to 
equation ( 1 ). 

The pressure distributions found by adding equations ( 1 ) and (2) are com- 
pared in figure 3 with experimental results from reference 8 . For M e = 1.322 
and Re * 9.4 x 105, other parameters are calculated as e * .247, u T = .050, 

6 = .0182, and c^ * .0032. The shock-wave position found from the data is the 
origin for the pressure at 15 mm., where y « 4.1. The origin for the wall 
pressure is shifted upstream through Ax = 2.87, equal to the displacement of the 
shock wave calculated using the fact that the shock-wave slope is proportional 
to the y-component of velocity at x = 0. In the experiments the wall was a 
plate having longitudinal curvature which can be inferred from the measured 
pressures immediately behind the shock wave. The change AP/P 0 in a y-distance 
of 15 mm. was taken as .015, giving K * .21. The theoretical curves also in- 
clude a constant streamwlse pressure gradient caused by cross-section area 
change, estimated from the data as AP/P 0 = .03 in a distance of 50 mm. The 
parameter II was taken equal to the constant-pressure value II = 0.5, which leads 
to a predicted location of the undisturbed sonic line at y = .36. If instead 
II = 0.4, the prediction is y = .30, whereas the measured value is y = .29; this 
change in II would influence the pressures only slightly. The wall pressures are 
in excellent agreement for x greater than about 3 ; the region x < 3 requiring 
solution of nonlinear equations is relatively large in this case, because the 
Reynolds number is not high enough for the upstream sor.c line to be really 
close to the wall . The very beginning of the rise in wall pressure can be ex- 
pressed by an exponential function with known rate of decay but unknown multipli- 
cative constant, as shown in figure 3 with origin chosen tentatively for good 
agreement with the data. By comparison with umerical results of ref. 2, it 
should be possible to estimate this constant. The pressures at 15 mm. from the 
wall ar in good agreement except for small x, where the error probably arises 
from approximating the shock-wave location by x = 0 ; adding the appropriate 
second-order term is expected to improve the agreement. 


WALL SHEAR 


The above description of mean-flow perturbations in terms of disturbances 
in an inviscid rotational flow does not remain valid as y -+■ 0. For y ■ Y /6 
“ °(u T e 1 ' 2 ) , it is easily found that turbulent stresses as well as pressure and 
inertia terms must be retained in the boundary-layer momentum equation. This 
thinner layer (still much thicker than the viscous wall layer) has been calle 
Reynolds stress sublayer in reference 1 and a blending layer in reference 2. 
Double expansions as u T -*■ 0 and e -*• 0 are assumed for the velocity components, 
pressure, and Reynolds stress, with the latter represented by a simple mixing- 
length model, chosen for analytical convenience Ln the belief that the results 
may not depend strongly on the model used. The wall shear stress has the form 

T w (x) “ 1 + a£ 4 -|a(a-l )£ 2 +• •• + u t t 1 (x) +eu T { (In f 1 / 2 u t ) '^(x) + t 11 (x)} 4 •• • ( 3 ) 

where a ■ constant. This expansion, like that for the pressure, requires modifi- 
cation when x is small. The y-component of velocity and the related displacement 


2/4 



effect of this sublayer are small enough that the previous calculation of pres- 
sure is unaffected, and the boundary condition v = 0 at y = 0 used for y = 0(1) 
is thereby verified With use of the solution for pressure evaluated as y - 0 , 
it is found that T^x) is proportional to the pressure perturbation, and there- 
fore decreases monotonically with increasing x; also t is constant. The term 
11 W is obtained as a lengthy expression including terms proportional to the 
pressure perturbation and a positive tenn proportional to In x (ref. 9). Since 

J/25 r L SS v Ur L P " tUrba , tl0n decreases to zero a s x increases (i.e., is small for 
m D * !)» it follows that t w has a minimum. Thus for typical values of 

M e and Re figure 4 shows that t w first decreases sharply and then rises again 

u A 5 3 ? iV6n *’ Tw decreases as E Increases and increases as Re increases. 
This behavior is in general agreement with existing experimental results. 

If it is assumed that equation (3) remains a good approximation even when 
changes in T w are not small as originally assumed, a criterion for incipient 

C ! n „ be p ™P° sed * For sample, curve III in figure 4 is drawn for 
Ke - 10 and M e = 1.149, corresponding to e = 0.12; if e is increased to 0.2, 
the minimum value of x w decreases to zero, and this is taken as an indication of 
incipient separation. Values of M e corresponding to this condition, namely that 
% - 0 and dx w /dx * 0 simultaneously, are plotted against Re* in figure 5, and 
show a slow increase in the shock strength required for separation as Re in- 
creases. As a first step toward studying flows with separation, a numerical 
solution was also carried out for a momentum equation appropriate for describing 
a slender low-speed separation bubble having length Ax * 0(1). Terms v/ere 
retained to represent the laminar and turbulent stresses and the pressure gra- 
dient, which was considered known from the previous solutions because the dis- 
placement effect of the bubble was sufficiently small. Velocity profiles 
obtained are plotted against y in figure 6. 



CONCLUDING REMARKS 


The favorable comparison with experiment shown in the figures for the latter 
part of the pressure rise seems to imply that the asymptotic representation 
given here, with the addition of second-order terms as noted, does correctly 
include the most important features of the local flow field. A study of other 
second-order terms is continuing, and an analytical --olution has been derived 
for the correction necessary in the flow behind as).'; wave in a circular pipe- 
a comparison with existing experimental results is .ng carried out. It Is 
anticipated that a comparison with numerical resu A of Melnick and Grossman 
wiil suggest a simple curve fit, in terms of a coou.inate nondimensional with 
the distance to the sonic line, to fill in the * -.11 pressure in the region near 
the origin. It should now be possible to incorporate these analytical representa- 
tions for the local preset re and velocity into poter.tial-f low calculations of 
the transonic flow past an airfoil. The predictions for wall shear, and for the 
Mach number at which separation first occurs, show the expected trends, and 
sh'-jld be checked for quantitative accuracy against existing experimental data. 

It would seem especially important to continue the attempt at developing a 
rctional theory of separation for turbulent boundary layers. 


I 
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Figure 5.- Predicted variation of Mach number with Reynolds nunber at 

incipient separation. 
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